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Abstract 

Fuel  cells  with  a  polymer  electrolyte  membrane  have  been  receiving  more  and  more  attention.  Modeling  plays  an  important  role  in  the 
development  of  fuel  cells.  In  this  paper,  the  state-of-the-art  regarding  modeling  of  fuel  cells  with  a  polymer  electrolyte  membrane  is  reviewed. 
Modeling  has  allowed  detailed  studies  concerning  the  development  of  these  cells,  e.g.  in  discussing  the  electrocatalysis  of  the  reactions  and 
the  design  of  water-management  schemes  to  cope  with  membrane  dehydration.  Two-dimensional  models  have  been  used  to  represent  reality, 
but  three-dimensional  models  can  cope  with  some  important  additional  aspects.  Consideration  of  two-phase  transport  in  the  air  cathode  of  a 
proton  exchange  membrane  fuel  cell  seems  to  be  very  appropriate. 

Most  fuel  cells  use  hydrogen  as  a  fuel.  Besides  safety  concerns,  there  are  problems  associated  with  production,  storage  and  distribution  of 
this  fuel.  Methanol,  as  a  liquid  fuel,  can  be  the  solution  to  these  problems  and  direct  methanol  fuel  cells  (DMFCs)  are  attractive  for  several 
applications.  Mass  transport  is  a  factor  that  may  limit  the  performance  of  the  cell.  Adsorption  steps  may  be  coupled  to  Tafel  kinetics  to  describe 
methanol  oxidation  and  methanol  crossover  must  also  be  taken  into  account.  Extending  the  two-phase  approach  to  the  DMFC  modeling  is  a 
recent,  important  point. 
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1.  Introduction 

Modeling  plays  an  important  role  in  the  development  of 
fuel  cells,  because  it  allows  a  better  comprehension  of  the 
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parameters  affecting  the  performance  of  single  fuel  cells  and 
fuel  cell  systems  [1]. 

Fuel  cells  with  a  polymer  electrolyte  membrane  have 
been  receiving  more  and  more  attention.  The  main  char¬ 
acteristics  of  this  kind  of  fuel  cell  are  the  operation  with¬ 
out  the  generation  of  pollutants,  less  corrosion  problems, 
high  power  density  and  low  temperature  start-up.  However, 
there  are  several  problems  to  be  solved,  such  as  mem- 
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Nomenclature 

a 

effective  catalyst  area  per  volume  (cm-1) 

^ht 

heat  transfer  area  per  unit  length  (cm) 

b 

Tafel  slope  (V) 

^cell 

parameter  (V  dec-1) 

Ca 

dimensionless  capacitance  for  the  agglomerate 
network 

^ch3oh 

CF 

methanol  feed  concentration  (mol  dm-3) 

Cf 

concentration  of  membrane  fixed-charge- site 
species  (mol  cm-3) 

Ck 

molar  concentration  of  k  (mol  cm-3) 

cw 

concentration  of  water  in  the  membrane 
(mol  cm-3) 

Cdi 

double  layer  capacitance  (F) 

Ci 

parameter  (V) 

c2 

parameter  (cm2  A-1) 

d 

channel  height  (cm) 

Df 

effective  mass  diffusivity  (cm2  s-1) 

Dx 

corrected  diffusion  coefficient  (cm2  s-1) 

Dij 

diffusion  coefficient  (cm2  s-1) 

E° 

parameter  (V) 

E* 

parameter  (V) 

f 

frequency  (rpm) 

F 

molar  flow  rate  of  methanol  (mol  s-1) 

F'j' 

total  molar  flow  rate  of  gas  (molm-1  s-1) 

h 

channel  width  (cm) 

h% 

height  of  the  active  catalyst  layer  (m) 

hf 

height  of  the  flow  channel  (m) 

hp 

height  of  the  porous  backing  (m) 

ht 

heat  transfer  coefficient  (W  m-2  K-1) 

Hem 

location  of  cathode  catalyst  layer  (cm) 

AH 

enthalpy  gradiente  (J  mol-1) 

i 

current  density  (A  cm-2) 

id 

limiting  current  density  (A  cm-2) 

io 

exchange  current  density  (A  cm-2) 

I 

current  density  (A  cm-2) 

h 

current  density  at  the  limiting  current  density 
(A  cm-2) 

2 

limiting  current  density  due  to  limiting  oxygen 
diffusion  (A  cm-2) 

4 

parasitic  current  density  at  cathode  (A  cm-2) 

j 

(-d1/2 

jMeOH 

methanol  mass  flux  (kg  cm  2  s  ]) 

jo 

exchange  current  density  referenced  to  pure 
oxygen  at  1  atm  (A  cm-2  atm-1) 

J 

current  density  (A  cm-2) 

kp 

hydraulic  permeability  (cm2) 

k<p 

electrokinetic  permeability  (cm2) 

F 

cell  channel  length  (cm) 

mu+ 

molality  of  H+  (molg-1) 

Mm 

equivalent  weigh  of  membrane 

Mi 

molar  flow  rate  of  i  (mol  s-1) 

^MeOH 

molecular  weight  of  methanol  (kg  mol-1) 

d k/dx  water  condensation  or  evaporation 
(mol  cm-1  s-1) 

^drag  electro-osmotic  drag  coefficient 

Niyfi  molar  flux  of  i  in  k  channel  (mol  cm-2  s-1) 

N(  molar  flux  of  i  (mol  cm-2  s-1) 

Nw  water  molar  flux  (mol  cm-2  s- 1 ) 

AW;C  water  molar  flux  produced  at  cathode 
(mol  cm-2  s-1 ) 

p  hydraulic  pressure  (atm) 

P  pressure  (atm  or  Pa) 

Pc  critical  pressure  (atm) 

Pw  partial  pressure  of  water  in  the  membrane  (atm) 
R  resistance  (£2  cm2) 

R\f  low  frequency  resistance  (£2) 

Rm  membrane  resistance  (£2  cm2) 

Rs  solution  resistance  (£2) 

Pi  dimensionless  resistance  in  the  agglomerate 
network 

.v  liquid  water  saturation 

SH  cell  height  (cm) 

SW  cell  width  (cm) 

S  flooding  parameter 

Su  momentum  source  term  (kg  cm-3  s-2) 

Sk  species  source  term  (mol  cm-3  s-1) 

Sep  potential  source  term  (S  V  cm-3) 

tm  membrane  thickness  (cm) 

t\  anode  thickness  (cm) 

tc  cathode  thickness  (cm) 

T  temperature  (°C  or  K) 

Pa  temperature  of  the  anode  stream  (°C) 

Tc  temperature  of  the  cathode  stream  (°C) 

Ps  temperature  of  the  solid  phase  (°C) 

u  velocity  vector  (cm  s  - 1 ) 

f/ht  overall  heat  transfer  coefficient 

(J  s-1  cm-2  °C-1) 

Us  (1  +  S/M)  (where  S/M  is  the  molar  ratio  of 
steam  to  methanol) 
v  velocity  (cm  s-1) 

Veen,  f/cell,  E  cell  potential  (V) 

Voc  open  cell  potential  (V) 

W  weight  of  catalyst  (kg) 

v  distance  (cm  or  m) 

Xi  mole  fraction  of  i 

v0,c  mole  fraction  of  oxygen  at  cathode 
viiq  mole  fraction  of  liquid  water 

vmet  conversion  ratio  of  methanol 

y  distance  (cm  or  m) 

Ys  C/(n(0)  +  r2(0))  (molkg-1  s_1) 

z  distance  (cm  or  m) 

zH+  charge  number  of  H+ 

Zf  charge  number  of  membrane  fixed-charge- site 
species 

Z  impedance  in  the  electrolyte  (£2) 


s 

SH 

SW 

s 

Su 

S* 

s0 

tm 

tA 

tc 

T 

T* 

Tc 

Ts 

u 

Uht 


Ao,c 

-^liq 

-^met 

J 

Ys 

z 

Zr+ 

Zf 
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^electrode  total  impedance  (Q) 

Zs  (Uri(0)  -  2r2(0))  (molkg-1  s_1) 

Greek  letters 

o' a  anodic  transfer  coefficient 

ctfc  cathodic  transfer  coefficient 

s  porosity 

£^+  cathode  gas  porosity 

0  potential  (V) 

0  dc  potential-independent  agglomerate 

diffusion  parameter 
r  thin  film  diffusion  parameter 

0  overpotential  (V) 

X  water  content  in  the  membrane 

ji  viscosity  (gem-1  s-1) 

fjc  empirical  constant 

/xeff  effective  viscosity  (kg  m- 1  s~ 1 ) 

/xH+  chemical  potential  of  H+  (J  mol  - 1 ) 

liv  reduced  overpotential 

p  density  (kg  cm-3) 

Pdry  density  (g  cm-3) 

creff  effective  ionic  conductivity  (S  cm- 1 ) 

crm  membrane  conductivity  ( Q  ~ 1  cm- 1 ) 

m  angular  frequency  (rad  s  - 1 ) 

z&a  characteristic  agglomerate  frequency  (rad  s  - 1 ) 

zuf  characteristic  thin  film  frequency  (rad  s  - 1 ) 

0  ac  potential-dependent  agglomerate  diffusion 

parameter 

f  stoichiometric  flow  ratio  (amount  of  reactant 

in  the  chamber  feed  divided  by  the  amount 
required  by  the  electrochemical  reaction) 
based  on  the  reference  current  density  of 
1  A  cm-2. 


brane  dehydration,  because  of  the  water  transport  in  the  fuel 
cell,  the  high  catalyst  cost  [2]  and  the  poisonous  effect  of 
contaminants. 

The  optimal  modeling  approach  differs  for  each  applica¬ 
tion.  It  is  important  to  define,  initially,  the  area  of  interest 
of  the  model.  It  can  be  at  the  fundamental  cell  level,  includ¬ 
ing  the  inlet  channels,  the  electrodes  and  the  membrane,  on 
a  higher  level  where  individual  fuel  cells  are  gathered  in  a 
stack,  or  on  the  complete  system  level  consisting  of  a  fuel 
cell  stack  with  its  compressor,  pumps,  humidifiers,  etc.  [1]. 
In  Haraldsson  and  Wipke  [1],  the  analysis  is  directed  to  the 
evaluation  of  fuel  cell  system  models.  In  this  review,  we  will 
mainly  be  concerned  with  modeling  at  the  fundamental  cell 
level. 

A  mechanistic  fuel  cell  model  is  based  on  transport  phe¬ 
nomena  and  electrochemical  relationships.  On  the  other  hand, 
there  are  empirical  models,  developed  specifically  for  each 
application  and  set  of  operating  conditions  (such  as  in  Kim  et 
al.  [3]  and  Squadrito  et  al.  [4]).  There  are  still,  semi-empirical 


models  in  which  the  introduction  of  some  phenomenologi¬ 
cal  knowledge  into  empirical  relations  can  improve  the  pre¬ 
dictability  of  the  empirical  model.  When  a  tool  for  detailed 
studies  is  desired,  as  many  as  possible  mechanistic  theoretical 
considerations  must  be  taken  into  account. 

The  operating  regime  associated  with  the  model  can  be 
either  steady-state  or  transient.  Defining  spatial  dimension 
and  complexity  is  also  important.  Validation,  in  turn,  is  impor¬ 
tant  to  guarantee  that  the  model  is  a  useful  and  reliable  tool. 


2.  Fuel  cell  modeling 

Springer  and  Raistrick  [5]  took  into  account,  in  the  model¬ 
ing  of  a  porous  electrode  with  oxygen  diffusion,  the  “flooded- 
agglomerate”  approach  of  Giner  and  Hunter  [6]  coupled  to  a 
thin-film  diffusion  process.  The  electrode  structure  was  con¬ 
sidered  to  be  a  macro-microporous  interconnected  system, 
with  hydrophobic  regions  allowing  reactant  gas  access  to  the 
surface  of  agglomerate  regions.  The  agglomerate  consists  of 
regions  of  catalyst  containing  carbon  and  electrolyte,  and  is 
separated  from  the  hydrophobic  gas  regions  by  a  thin  film 
of  electrolyte.  Reactant  gas  can  dissolve  in  the  electrolyte 
and  diffuse  to  the  catalyst  sites.  Transfer  functions  developed 
allow  a  combination  of  a  simplified  agglomerate  representa¬ 
tion  with  the  thin  film  element,  facilitating  numerical  fitting 
of  the  model  to  measured  impedance  data.  Perez  et  al.  [7]  used 
Springer  and  Raistrick  approach  (for  simulation  and  analysis 
of  dc  and  ac  polarization  results)  in  studies  on  the  electrocatal¬ 
ysis  of  the  oxygen  reduction  reaction  on  platinum  on  carbon 
thin  porous  coating  rotating  disk  electrode  (TPC/RDE).  The 
TPC/RDE  is  an  experimental  apparatus  conceptually  similar 
to  a  flooded-agglomerate,  given  by  the  electrode  structure, 
in  contact  with  a  thin  film  of  thickness  determined  by  the 
rotation  rate  of  the  electrode.  Steady-state  polarization  (fit¬ 
ted  using  Eq.  (la))  and  the  impedance  (fitted  using  Eq.  (lb)) 
results  showed  a  duplication  of  the  Tafel  slope  (Fig.  la)  in 
alkaline  medium  (in  acid  medium,  there  is  a  duplication  first 
and  then  a  quadruplication,  Fig.  lb).  The  duplication  in  alka¬ 
line  medium  and  the  first  duplication  in  acid  medium  are  due 
to  structural  effects  in  the  catalyst  layer  of  the  electrode  [7]: 

i  exp/xr(tanh[0  exp(/xr/2)]  )/[<£>  exp(/zr/2)] 

—  = -  (la) 

io  1  +  Eexp/xr(tanh[0exp(/xr/2)])/[0exp(/xr/2)] 


Z  =  {exp  /zr(/0 /£>)/(  1  +  r  exp/xr/R/) 

x[l/(/?i  +(l/((l/(((l/2)(sech20  exp(Mr/2)) 
+(l//?l))_1  -  R\))  +  jca(m/ma)))) 

-((T  expMr/fli)(l  +  ( <P 4  exp(2Mr)/16))1/4 
x(tanh  tf'/<f02(tanh  V jUn/m)/  flG/m)))/ 

(1  +  r  exp ixr(^I/2/(02  exp  /xr))(tanh  'I'/'I') 

(tanh  J j(uj/ujf)/  J j{zu/mi)))]}~1  (lb) 
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Fig.  1.  (a)  Mass-transport  corrected  Tafel  plots  for  oxygen  reduction  on 
40%  Pt/C  in  1.0  M  NaOH.  (x)  log(l/7?if)  vs.  E;  (O )  log[(/d  x  1)1  (Id  -  /)]  vs. 
E\  (full  line)  fitted  plot  using  Eq.  (la)  with  r  — >  0./=2500rpm  (reprinted 
from  [7],  Copyright  1998,  with  permission  from  Elsevier  Science),  (b)  Mass- 
transport  corrected  Tafel  plots  for  oxygen  reduction  on  40%  Pt/C  in  0.5  M 
H2SO4.  (x)  log(l/flif)  vs.  E\  (O )  log[(/d  x  /)/(/d  —  /)]  vs.  E\  (full  line)  fitted 
plot  using  Eq.  (la)  with  T  — >  0./=  2500  rpm  (reprinted  from  [7],  Copyright 
1998,  with  permission  from  Elsevier  Science). 


In  Springer  et  al.  [8],  the  distribution  of  water  in  a  polymer 
electrolyte  fuel  cell  is  calculated  by  considering  the  water 
flow  through  two  inlet  channels,  two  gas-diffusion  electrodes 
and  a  Nation  membrane  (Fig.  2). 

Water  enters  the  fuel  cell  assembly  as  a  component  of 
the  humidified  fuel  and  oxidant  streams,  and  by  generation 
associated  with  the  cell  reaction.  Diffusion  of  water  vapor  and 
hydrogen  occurs  through  the  anode,  while  ternary  diffusion  of 
water  vapor,  oxygen  and  nitrogen  occurs  through  the  cathode. 
These  diffusion  processes  through  the  porous  electrodes  are 
calculated  starting  from  the  Stefan-Maxwell  equation  (Eq. 
(2)): 


d  Xj 

d  z 
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Fig.  2.  Schematic  diagram  of  fuel  cell  model  (reprinted  from  [8],  Copyright 
1991,  reproduced  by  permission  of  The  Electrochemical  Society,  Inc.). 


Equilibrium  of  water  is  assumed  at  the  electrode/memb¬ 
rane  interface.  The  net  flux  of  water  through  the  membrane 
considers  the  electro-osmotic  and  the  diffusion  driving  forces 
for  water  (/Vw  —  ^drag(2/VW5C)(AV22)  (pdry//\Tm)Z)^(dA/dz)). 

When  the  material  balances  are  solved,  the  membrane  resis¬ 
tance  Rm  =  Jjm(l /cr(A))  d z,  the  O2  concentration  at  the  cat¬ 
alyst  interface  and  the  cathode  overpotential  (by  using  the 
Tafel  expression  J=joPc(x0?c/(l  —  jcnq))  Qxp[0.5Fr]/(RgT)]) 
are  determined.  So,  a  potential  versus  current  curve  for 
the  fuel  cell  can  be  determined  (VCeii  =  VoC  —  0  —  JRm).  The 
active  catalyst  layer  is  assumed  to  exist  as  a  thin  plane  at 
the  electrode/membrane  interface  (as  it  is  considered  in  the 
majority  of  the  complete  cell  models);  the  rest  of  the  elec¬ 
trode  serves  only  as  a  gas  diffusion  region.  The  Springer  et 
al.  model  predicted  the  increase  in  membrane  resistance  with 
increased  current  density  and  pointed  out  the  advantage  of 
thinner  membranes  in  reducing  this  problem. 

In  Bernardi  and  Verbrugge  [9],  a  similar  modeling  of 
a  polymer  electrolyte  fuel  cell  is  presented.  The  cell  con¬ 
sists  of  a  membrane  sandwiched  between  two  gas-diffusion 
electrodes  (membrane  and  electrodes  assembly,  MEA),  hot 
pressed  and  placed  between  two  current  collectors.  Humidi¬ 
fied  gaseous  H2  enters  the  anode  gas  chamber,  is  transported 
through  the  porous  gas  diffusion  layer  and  dissolves  in  the 
electrolyte  phase  at  the  anode  catalyst  layer  where  it  is  oxi¬ 
dized.  The  gaseous  reactant  O2,  mixed  with  nitrogen  and 
water  vapor,  enters  into  the  cathode  gas  chamber,  is  trans¬ 
ported  through  the  porous  gas  diffusion  layer  and  dissolves 
in  the  electrolyte  phase  at  the  cathode  catalyst  layer.  Pro¬ 
tons  transported  through  the  membrane  participate  in  the 
reduction  of  the  dissolved  O2  at  catalyst  sites.  Bernardi  and 
Verbrugge  assume  that  the  system  operates  at  constant  tern- 
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perature.  It  is  also  assumed  steady-state  operation.  They 
consider  a  fully  hydrated  membrane  and  wet  pores  in  the 
gas  diffusion  layer  (assuming  no  interactions  between  liquid 
and  gas  flows).  The  equations  constituting  the  mathematical 
model  of  the  cell  are  derived  from  basic  phenomenological 
relationships  (such  as  a  form  of  Schlogl’s  velocity  equation 
(Eq.  (3a)),  the  Butler- Volmer  equation  for  charge  transfer 
(Eq.  (3b))  and  the  Stefan-Maxwell  equation  for  gas-phase 
transport): 


0 


These,  along  with  appropriate  conservation  principles  of 
physics  allow  deriving  the  coupled  differential  equations  for 
the  cell  model.  The  necessary  parameters  and  properties  were 
appropriately  considered  and  estimated.  For  the  base-case, 
results  do  not  show  cell  polarization  resulting  from  limitation 
on  the  transport  of  O2  (Fig.  3a).  However,  by  decreasing  the 
cathode  gas  porosity,  the  limited  transport  of  oxygen  through 
the  cathode  diffusion  layer  to  the  reaction  sites  may  increase 
the  cathode  overpotential  (Fig.  3b). 

It  also  showed  that  membrane  dehydration  could  cause 
limitations  on  the  operating  current  density.  In  this  case,  effi¬ 
cient  water-management  schemes  are  necessary. 

In  1993,  Fuller  and  Newman  published  a  paper  [10]  in 
which  it  is  considered  a  fuel  cell  operating  in  a  steady- 
state  regime  on  air  and  reformed  methanol.  The  transport 
of  species  in  the  polymer  electrolyte  was  described  by  a 
fundamental  equation  that  may  be  put  in  the  form  of  the 
Stefan-Maxwell  equation  and  combined  with  the  appropri¬ 
ate  material  balance  of  each  species.  The  Stefan-Maxwell 
equations  could  also  be  used  to  describe  multicompo¬ 
nent  diffusion  of  gases.  Electron-transfer  reactions  were 
described  by  Butler- Volmer  kinetics.  The  energy  released 
by  the  overall  chemical  reaction  was  split  into  electrical 
work,  change  in  the  enthalpy  (and  hence  in  the  tempera¬ 
ture)  of  the  flowing  gas  streams  and  external  heat  transfer 
(Eq.  (4)).  Heat  removal  is  also  important  in  the  fuel  cell 
operation: 

f  ht(T  -  ^Ambient)  dz  +  f  W  dz  +  FT  AH  =  0  (4) 

Jo  Jo 

In  1993,  Nguyen  and  White  were  also  concerned  with  the 
water  and  heat  management  in  proton  exchange  membrane 
fuel  cells,  which  are  essential  features  for  obtaining  high 
power  densities.  A  water  and  heat  management  model  was 
developed  [11]  and  used  to  investigate  the  effectiveness  of 
different  humidification  designs.  The  model  is  a  steady-state, 
two-dimensional  one.  The  model  regions  (Fig.  4)  consist  of 
flow  channels  on  both  sides  of  the  membrane/electrode  (one 
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Fig.  3.  (a)  Model  calculations  of  the  contributions  to  fuel  cell  potential  losses 
as  a  function  of  operating  current  density  for  base  case  conditions  (reprinted 
from  [9],  Copyright  1992,  reproduced  by  permission  of  The  Electrochem¬ 
ical  Society,  Inc.),  (b)  Model  calculations  of  the  fuel  cell  potential  as  a 
function  of  operating  current  density  for  two  values  of  cathode  gas  porosity 
gg+  (reprinted  from  [9],  Copyright  1992,  reproduced  by  permission  of  The 
Electrochemical  Society,  Inc.). 


system  for  the  oxygen  cathode  and  another  one  for  the  hydro¬ 
gen  anode). 

The  model  takes  into  account  mass  transport  of  water 
and  gaseous  reactants  across  the  membrane/electrode  and 
along  the  flow  channels  and  heat  transfer  between  the  solid 
phases  and  the  gases  along  the  flow  channels.  The  change 
in  the  number  of  moles  (quantified  by  the  material  balances) 
of  a  single-phase  species  along  the  channel  length  is  due 
to  the  normal  flux  into  or  out  of  the  membrane/electrode 
(6Mi/dx  =  hNitytk).  For  water,  the  material  balance  also  takes 
into  account  the  possibility  of  condensation  or  evaporation 
(quantified  by  the  term  (d k/dx)).  It  is  assumed  that  water 
leaves  the  channel  and  enters  the  electrode  in  the  form  of 
vapor  only.  Fiquid  water  is  assumed  to  exist  in  the  form 
of  small  droplets  with  negligible  volume.  The  electrodes 
are  considered  ultrathin  (gas  diffusion  through  the  electrode 
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Fig.  4.  Schematic  of  modeled  regions  (reprinted  from  [11],  Copyright  1993, 
reproduced  by  permission  of  The  Electrochemical  Society,  Inc.). 

porous  layer  is  neglected).  The  molar  flux  of  each  compo¬ 
nent  is  proportional  (depending  on  the  stoichiometry  of  the 
electrochemical  reactions)  to  the  local  current  density,  which 
varies  along  the  channel  length  as  the  membrane  conduc¬ 
tivity  and  the  overpotential  of  the  electrodes  change.  The 
net  water  flow  in  the  membrane  is  associated  to  migration 
(water  molecules  carried  by  protons)  and  diffusion  effects. 
Adequate  expressions  to  estimate  important  parameters,  such 
as  the  water  diffusion  coefficient  and  the  electro-osmotic  drag 
coefficient,  were  used.  The  energy  balance  equation  for  the 
anode  and  cathode  gaseous  streams  takes  into  account  water 
condensation  and  evaporation  (proportional  to  d Mlw  k/dx)  in 
the  flow  channel  and  heat  transfer  due  to  the  temperature  gra¬ 
dient  between  the  solid  phase  and  the  interior  of  the  channels. 
The  cell  potential  is  calculated  accounting  for  the  membrane 
resistance  and  the  electrode  polarizations.  Modeling  results 
showed  that  at  high  current  densities,  the  back  diffusion  of 
water  from  the  cathode  side  of  the  membrane  is  insufficient  to 
keep  the  membrane  hydrated.  Water  injection  into  the  anode 
flow  channel  can  improve  the  cell  performance  for  supplying 
the  humidification  to  maintain  the  anode  hydrated,  and  also 
for  contributing  to  the  heat  removal.  When  air  is  used  (instead 
of  pure  oxygen),  the  cathode  stream  must  also  be  humidified. 

In  the  previously  cited  Nguyen  and  White’s  work  [11], 
the  temperature  of  the  solid  phase  (flow  channel  plates, 
electrodes  and  membrane)  was  considered  uniform  and  con¬ 
stant.  In  practice,  a  temperature  gradient  in  the  solid  layers 
along  the  flow  channels  of  the  fuel  cell  may  exist  (the  term 
£4t<2ht(^solid  —  Tchannel)  accounts  for  the  energy  transported 
from  the  fluid  to  the  solid  layer),  and  this  gradient  can  affect 
the  temperature  distribution  of  the  fluid  in  the  flow  channels 
and  also  the  hydration  state  of  the  membrane  along  the  flow 
path.  To  better  deal  with  these  problems,  a  new  model  was 


1  sw  i 

Fig.  5.  Schematic  diagram  of  modeling  region  (reprinted  from  [12],  Copy¬ 
right  1998,  reproduced  by  permission  of  The  Electrochemical  Society,  Inc.). 


developed  [12],  based  on  Nguyen  and  White  [1 1],  to  describe 
the  temperature  distribution  in  the  solid  phase  along  the  flow 
path  of  the  anode  and  cathode  faces  of  the  fuel  cell  (Fig.  5). 

The  model  was  extended  to  include  the  convective  water 
transport  across  the  membrane  by  a  pressure  gradient  (quan¬ 
tified  by  the  term  cw(kp//x)(dJPw/dy))  and  the  heat  removal 
by  heat  exchangers.  A  pressure  gradient  between  the  cath¬ 
ode  and  anode  can  increase  the  water  back  transport  across 
the  membrane  (although  it  can  be  demonstrated  that  the  fuel 
cell  pressurization  involves  an  energy  cost  that  hardly  com¬ 
pensates  the  gain  in  performance — De  Souza  and  Gonzalez 

[13] ),  while  using  appropriate  heat  exchangers  can  be  amply 
effective  for  heat  removal  and  distribution.  Dannenberg  et  al. 

[14]  presented  a  modeling  work  according  to  an  approach 
very  similar  to  those  of  the  previously  cited  Nguyen  and  col¬ 
laborators’  works. 

In  the  work  of  Sena  et  al.  [15],  a  simplified  model  is  used 
to  describe  the  water  transport  in  a  polymer  electrolyte  mem¬ 
brane  of  a  H2-O2  fuel  cell  operating  at  low  temperatures; 
a  diffusion  mechanism  is  considered  for  water  back  trans¬ 
port.  Combined  with  the  Tafel  equation,  it  led  to  an  analytical 
expression  for  the  description  of  the  cell  potential  as  a  func¬ 
tion  of  the  current  density.  Additionally,  limiting  effects  due 
to  oxygen  diffusion  were  included  (Eq.  (5)): 

£  =  £°  -  Mog/ +  Mog  (l  -  -M  (5) 

Results  showed  that  for  Nation  115  and  117  membranes, 
water  transport  through  the  membrane  is  a  limiting  factor.  For 
Nation  112,  on  the  other  hand,  oxygen  diffusion  effects  (in 
the  gas  diffusion  electrode)  dominate  the  potential-current 
relationship. 

In  Gurau  et  al.  [  1 6] ,  a  two-dimensional  model  for  the  entire 
sandwich  of  a  H2-O2  proton  exchange  membrane  fuel  cell 
was  developed.  Standard  momentum,  energy  transport,  con¬ 
tinuity  and  species  concentration  equations  were  solved  in 
the  gas  channels,  while  the  porous  media  model  was  used  for 
the  equations  describing  transport  phenomena  in  the  mem¬ 
brane  and  the  gas  diffusion  electrodes  (in  particular,  with 
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Fig.  6.  Effect  of  the  inlet  hydrogen  mole  fraction  on  cell  polarization  curves 
(reprinted  from  [17],  Copyright  2000,  reproduced  by  permission  of  The 
Electrochemical  Society,  Inc.). 


the  velocity  field  provided  by  a  generalized  form  of  Darcy’s 
law,  obtained  from  the  standard  Navier-Stokes  equations). 
By  solving  the  transport  equations,  as  well  as  relationships 
for  the  electrochemical  reactions,  polarization  curves  under 
various  operating  conditions  were  obtained.  In  contrast  to 
the  approach  of  Gurau  et  al.,  which  used  separate  differential 
equations  for  different  subregions,  Um  et  al.  [17]  considered 
a  single-domain  approach,  in  which  a  single  set  of  equa¬ 
tions  valid  for  all  subregions  is  used.  The  fuel  cell  operation 
considering  isothermal  conditions  is  described  by  conserva¬ 
tion  principles  of  mass,  momentum,  species  and  charge  (Eqs. 
(6a)-(6d)): 


9(gp)  _ 
dt 

d(spu) 

dt 

d(sck) 

dt 

0  =  V  . 


—  V  •  ( spu ) 

=  —  V  •  (spun)  —  sV P  +  V  •  (£/xeffVw)  +  Su 

=  -V  •  (, euck) )  +  V  •  (DfVck)  +  Sk 
(ffeffV0)  +  S0 


(6a) 

(6b) 

(6c) 

(6d) 


The  momentum  source  term  is  used  to  describe  Darcy’s 
drag  for  flow  through  the  porous  electrodes  and  the  mem¬ 
brane.  The  species  and  potential  source  terms,  in  their  turn, 
are  related  to  generation  (or  consumption)  of  species  and  the 
creation  of  electric  current  (kinetic  expressions  were  derived 
from  the  Butler- Volmer  equation).  The  boundary  conditions 
are  required  only  at  the  external  surfaces  of  the  spatial  domain 
due  to  the  single-domain  approach  used.  The  conservation 
equations  were  discretized  and  solved  by  using  a  general- 
purpose  CFD  (computational  fluid  dynamics)  algorithm.  The 
model  could  show  (Fig.  6)  that  in  the  presence  of  hydrogen 
dilution  in  the  fuel  stream,  hydrogen  is  depleted  at  the  reac¬ 
tion  surface,  resulting  in  substantial  anode  mass  transport 
polarization  (this  study  can  be  of  interest  when  considering 
reformed  gas  as  the  anode  feed). 


In  a  subsequent  work,  Um  and  Wang  [18]  present  a 
three-dimensional  model  (Fig.  7)  to  cope  with  some  addi¬ 
tional  aspects,  such  as  the  study  of  a  mass- transport  limi¬ 
tation  not  discussed  in  Um  et  al.  [17],  the  rib  effect.  The 
two-dimensional  simulation  overpredicts  the  current  den¬ 
sity,  because  the  landing  area  obstructing  the  gas  transport 
is  neglected  [18].  The  design  of  an  innovative  flow  field, 
called  interdigitated,  was  also  studied.  Yi  and  Nguyen  [19] 
attempted  to  develop  a  two-dimensional  model  for  this  flow 
field.  De  Souza  and  Gonzalez  [13]  demonstrated  experimen¬ 
tally  the  advantages  of  the  interdigitated  flow  field.  In  this 
flow  field,  the  reactant  gas  is  forced  by  convection  to  pass 
through  the  porous  gas  diffusion  layer  in  order  to  leave  the 
system  (inlet  and  outlet  channels  are  not  connected). 

This  forced  convection  mechanism  facilitates  the  trans¬ 
port  of  reactants  and  products.  Two-dimensional  models  are 
limited  in  their  capability  to  represent  the  three-dimensional 
nature  of  the  interdigitated  flow.  The  three-dimensional  sim¬ 
ulations  are  again  based  on  a  single-domain  fuel  cell  model. 
The  model  consists  of  non-linear,  coupled  partial  differential 
equations  representing  the  conservation  of  mass,  momentum, 
species  and  charge  with  electrochemical  reactions.  Again,  the 
conservation  equations  were  discretized  and  solved  by  using 
a  general-purpose  CFD  algorithm. 

In  Dutta  et  al.  [20],  two-phase  transport  is  also  neglected. 
In  Janssen  [21],  a  two-phase  water  transport  model  is  pre¬ 
sented,  but  it  results  in  equations  that  are  difficult  to  imple¬ 
ment  in  a  general  purpose  CFD  algorithm.  In  Wang  et  al. 
[22],  the  so-called  Multiphase  Mixture  Model  [23,24]  is 
applied  to  study  the  two-phase  water  transport  in  the  porous 
air  cathode  of  a  proton  exchange  membrane  fuel  cell.  Water 
condensation  in  the  porous  cathode  is  a  usual  phenomenon 
at  high  current  densities.  As  the  current  density  (and  the 
water  generation)  increases,  water  is  present  first  as  vapor 
only.  When  the  water  vapor  density  at  the  membrane/cathode 
interface  reaches  the  saturation  value  corresponding  to  the 
operating  temperature,  liquid  water  begins  to  appear.  With 
further  increase  in  the  current  density,  the  two-phase  zone 
expands  and  its  front  propagates  towards  the  flow  channel 
(Fig.  8). 

In  the  mixture  approach  [23,24],  the  differential  balances 
are  applied  to  a  phase  mixture  (and  appropriate  additional 
algebraic  relationships  allow  to  obtain  information  on  each 
single  phase).  The  momentum  balance  equation  reduces  to 
an  extended  Darcy’s  law  for  two  phase  flow  in  porous  media. 
The  equation  for  the  transport  of  species  is  applicable  to 
oxygen  as  well  as  water.  The  mixture  model  is  mathemat¬ 
ically  equivalent  to  the  standard  two-phase  model  used  in  a 
large  number  of  applications  of  transport  in  porous  media, 
such  as  petroleum  exploitation  problems;  however,  the  mix¬ 
ture  model  is  easier  to  implement  [23,24].  In  Wang  et  al. 
[22],  the  cell  temperature  was  considered  constant,  and  the 
model  isothermal  (however,  an  energy  equation  may  improve 
the  analysis — Berning  et  al.  [25]).  Inside  the  two-phase 
zone,  the  transport  of  liquid  water  is  controlled  by  capillary 
action. 
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Fig.  7.  Three-dimensional  schematic  diagram  (reprinted  from  [18],  Copyright  2004,  with  permission  from  Elsevier  Science).  Straight  and  interdigitated  flows 
are  studied  by  changing  the  flow  inlet  and  outlet  locations. 


In  Pisani  et  al.  [26],  a  semi-empirical  model  is  presented, 
which  predicts  the  cell-potential  versus  current  density  (Eq. 
(7)).  The  derivation  had  the  objective  of  characterizing 
the  largest  number  of  mechanism-based  fitting  coefficients, 
instead  of  purely  empirical  ones  (such  as  in  Kim  et  al.  [3]  and 
Squadrito  et  al.  [4]).  It  was  based  on  the  fact  that  the  main 
non-linear  contributions  to  the  cell  potential  drop  in  H2-air 


Fig.  8 .  Liquid  water  saturation  contours  in  the  porous  cathode  and  flow  chan¬ 
nel  at  the  current  density  of  1.4  A  cm-2.  The  evaporation  front  separating 
the  two-phase  zone  from  the  single-phase  region  is  approximately  repre¬ 
sented  by  the  contour  of  s  =  0.01  (reprinted  from  [22],  Copyright  2001,  with 
permission  from  Elsevier  Science),  y  <0.07  cm,  flow  channel;  y  >0.07  cm, 
porous  cathode. 


fuel  cells  arise  in  the  cathode: 

Vcell  —  E°  —  RI  —  b  In  (/) 

+a  In  (l  -  (2)  s-Vc(i-(J/h))\  (7) 

In  Mazumder  and  Cole  [27],  a  three-dimensional  model  is 
presented,  which  carefully  extracts  the  good  elements  from 
the  existing  models  in  the  area.  In  some  aspects,  changes 
and  improvements  were  made.  The  modeling  is  based  on 
the  solution  of  conservation  equations.  Electrochemical  reac¬ 
tions  inside  the  porous  electrodes  are  described  by  using 
Butler- Volmer  equations.  The  catalyst  layer  is  treated  as  a 
finite  sized  region  rather  than  a  thin  plane  or  a  boundary.  The 
resistance  to  the  diffusion  of  the  gaseous  species  imposed 
by  the  electrolyte  (at  the  catalyst)  is  quantified  by  using  an 
effective  diffusion  coefficient  and  considering  that  the  overall 
porosity  of  the  active  catalyst  layer  is  smaller  than  that  of  the 
gas  diffusion  layer.  The  model  is  non-isothermal.  There  are 
some  fundamental  points  in  the  modeling. 

Porous  media  transport 

Conservation  equations  of  momentum  and  mass  in  a 
porous  media:  The  momentum  equation  is  of  the  type  of 
a  Darcy’s  equation. 

Conservation  equation  for  energy:  Electrical  work,  Joule 
heating  and  energy  interactions  because  of  phase  changes 
are  considered,  besides  the  irreversible  reaction  losses  due 
to  the  conversion  of  chemical  energy  into  heat. 

Species  conservation:  The  transport  of  species  in  gaseous 
regions  is  quantified  by  using  the  Stefan-Maxwell  equa¬ 
tions. 
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Polarization  Curve  for  323K  and  1/1  atm 


Fig.  9.  Comparison  of  experimentally  measured  performance  [29,30]  with 
numerical  predictions  with  and  without  liquid  water  formation  and  transport 
(reprinted  from  [28],  Copyright  2003,  reproduced  by  permission  of  The 
Electrochemical  Society,  Inc.). 

Electrochemical  reactions  and  coupling  to  mass  and  current 
transport 

Electrochemical  reactions  occur  at  the  interface  of  the 
catalyst  clusters  and  the  polymer  electrolyte.  It  is  where  the 
gas  comes  in  contact  with  the  catalyst  clusters  and  at  these 
points,  the  diffusion  flow  is  exactly  balanced  by  the  reaction 
flow.  The  charge  transfer  current  density  is  obtained  from 
the  Butler- Volmer  equation. 


term  to  account  for  membrane  ageing);  Maggio  et  al.  [35] 
(in  which  simulation  results  allow  a  comparison  of  different 
membranes);  Rowe  and  Li  [36]  (in  which  a  model  is  devel¬ 
oped  to  study  thermal  effects);  Ge  and  Yi  [37]  (in  which 
a  model  is  used  to  compare  counterflow  mode  to  coflow 
regime);  Yerramalla  et  al.  [38]  (in  which  a  fuel  cell  stack 
system  is  also  simulated). 

Most  of  the  fuel  cells  use  hydrogen  as  a  fuel.  Besides 
safety  concerns,  there  are  problems  associated  with  produc¬ 
tion,  storage  and  distribution  of  this  fuel.  Methanol,  as  a 
liquid  fuel,  can  be  the  solution  to  these  problems,  although  it 
presents  some  problems  of  toxicity  and  the  performance  of 
methanol  fuel  cells  is  lower  than  that  of  hydrogen  fuel  cells. 

Methanol  can  be  steam  reformed  to  CO2  and  H2;  at  the 
same  time,  some  of  the  methanol  decomposes  to  CO  and 
H2.  The  reforming  reaction  can  be  described  by  a  first  order 
kinetic  equation  r\  (in  relation  to  the  methanol  concentration) 
and  the  decomposition  one  by  a  zero  order  reaction  rate  r^ 
[39].  According  to  Rabou  [40],  the  problem  can  be  treated 
analytically  by  taking  the  methanol  conversion  ratio  as  the 
variable  of  interest,  and  expressing  the  first  order  kinetic 
equation  as  a  function  of  the  conversion  ratio  (through  the 
dependence  of  the  methanol  concentration  on  its  conversion 
ratio).  The  needed  catalyst  weight  for  a  given  methanol  flow 
and  conversion  fraction  can  be  calculated  as  a  function  of 
temperature  and  pressure  (Eq.  (8)): 


The  solution  of  the  set  of  conservation  equations  was  per¬ 
formed  within  the  framework  of  a  commercial  CFD  code. 

In  Mazumder  and  Cole  [28],  formation  and  transport  of 
liquid  water  in  the  proton  exchange  membrane  fuel  cell  is 
included.  The  model  is  based  on  the  multiphase  mixture 
model  (such  as  in  Wang  et  al.  [22]).  The  model  base  is 
an  additional  equation  for  the  formation  and  transport  of 
liquid  water.  Such  equation  includes  mainly  transport  associ¬ 
ated  with  electro-osmotic  drag,  influence  of  surface  tension, 
phase  change  (condensation  or  evaporation).  Liquid  water 
can  also  reduce  the  diffusivity  of  the  gaseous  species.  The 
“real”  porosity  for  the  gas  phase  species  is  called  wet  poros¬ 
ity  (and  is  a  function  of  the  liquid  saturation  of  the  pore,  which 
is  defined  as  the  volume  fraction  of  the  total  pore  space  occu¬ 
pied  by  the  liquid).  The  inclusion  of  liquid  water  formation 
and  transport  improved  the  predictive  capability  of  the  model 
(Fig.  9). 

Knudsen  diffusion  (which  occurs  when  the  pores  are  very 
tiny)  was  neglected,  and  it  is  worth  mentioning  that  Wohr  et 
al.  [31]  are  one  of  the  few  that  attempted  to  take  into  account 
Knudsen-diffusion  when  considering  the  diffusion  layer  of 
an  electrode,  but  no  discussion  is  presented. 

Additional  references  regarding  modeling  of  proton 
exchange  membrane  fuel  cells  are:  Eikerling  and  Kornyshev 
[32]  (in  which  a  homogeneous  model  for  the  cathode  is  con¬ 
sidered);  Kulikovsky  et  al.  [33]  (in  which  some  local  features 
generated  by  geometric  aspects,  e.g.  the  low  concentrations 
of  reactants  in  front  of  current  collectors,  are  studied);  Mann 
et  al.  [34]  (in  which  an  interesting  feature  of  the  model  is  a 


W  =  F 


In 


The  equation  can  also  be  inverted  to  calculate  the  methanol 
feed-flow  rate  necessary  for  a  determined  catalyst  weight  and 
a  given  conversion  fraction. 

Direct  methanol  fuel  cells  (DMFCs)  are  attractive  for 
several  applications.  The  structure  of  a  single  DMFC  is  an 
assembly  of  two  porous  electrodes  on  either  side  of  a  proton 
exchange  membrane.  Nowadays,  a  significant  performance 
has  been  achieved  by  using  an  anodic  catalyst  based  on  Pt- 
Ru,  while  for  the  oxygen  cathode  the  more  important  catalyst 
is  Pt. 

Jeng  and  Chen  [41]  present  a  mathematical  model  only 
for  the  anode  of  a  DMFC.  This  model  accounts  for  the  mass 
transport  in  the  flow  channel  (methanol  transport,  from  the 
feed  stream  to  the  diffusion  layer,  due  to  the  concentration 
gradient),  transport  in  the  diffusion  layer  (water  transport  and 
methanol  transport  by  diffusion  and  convection),  transport  in 
the  catalyst  layer  with  a  finite  thickness  (transport  similar  to 
that  in  the  diffusion  layer)  and  transport  in  the  membrane 
(water  transport,  by  diffusion  and  electro-osmotic  drag,  and 
methanol  transport  in  a  similar  way  as  in  the  diffusion  layer). 
The  electrochemical  reaction  rate  in  the  catalyst  layer  could 
be  described  by  using  a  kinetic  Tafel  expression.  In  Nord- 
lund  and  Lindbergh  [42],  an  agglomerate  model  (in  which 
each  agglomerate  is  modeled  as  a  sphere  with  a  content  of 
ionomer  and  catalyst  on  a  carbon  carrier)  is  used  and  kinetic 
equations  based  on  a  reaction  mechanism  (Eq.  (9))  are  derived 
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for  methanol  oxidation  in  the  anode  of  a  DMFC  (adsorption 
steps  coupled  with  Tafel  kinetics).  In  fact,  the  kinetic  model 
represents  a  simplified  mechanism  derived  from  a  more  com¬ 
plex  one  (Hamnett  [43]): 


CH3OH  &  CH3OHads 

*i 

CH3OHads  -X  COads  +  4H+  +  4e" 

H20  %  OHads  +  H+  +  e“ 

k> 

COads  +  OHads  -  CO2  +  1 1  +  e 

Since  heat  is  easily  transported  in  a  liquid  fuel,  tempera¬ 
ture  gradients  are  small  and,  therefore,  considered  negligible. 
Simulations  showed  that  mass-transport  limitations  within 
the  agglomerates  are  small  and  the  model  can  be  simpli¬ 
fied  (concentration  gradients  within  the  agglomerates  can  be 
neglected). 

In  Birgersson  et  al.  [44],  a  one-phase  model  for  mass, 
momentum  and  species  transport  in  the  anode  of  a  DMFC 
is  considered.  The  equations  and  boundary  conditions  were 
non-dimensionalized,  and  by  recognizing  that  the  ratio  of 
the  heights  of  the  flow  channel,  porous  backing  and  active 
catalyst  layer  to  the  anode  length  is  much  smaller  than  one 
(i.e.,  1,  /zp/L<jC  1  and  ha  <<C  1),  a  reduced  model  was 

derived.  The  reduced  model  requires  less  computing  time,  a 
factor  that  makes  it  suitable  for  inclusion  in  studies  where 
computational  time  may  be  a  central  matter.  Besides,  the 
study  of  the  effect  of  design  and  operating  parameters  on 
the  system  behavior  can  be  based  on  some  of  the  non- 
dimensional  parameters  that  arise  from  the  analysis. 

In  Scott  et  al.  [45],  a  DMFC  model  is  reported,  which 
accounts  for  the  influence  of  methanol  mass  transport  and  of 
the  CO2  flow  on  the  performance  of  the  fuel  cell.  The  model 
considers  the  hydrodynamics  in  the  flow  channels  associated 
with  the  production  of  CO2  and  methanol  mass  transport 
in  the  diffusion  electrode  and  the  polymer  membrane.  The 
model  also  incorporates  a  semi-empirical  equation  for  the 
open  circuit  potential  in  the  presence  of  methanol  crossover 
(in  practice,  a  model  of  the  effect  of  methanol  on  the  mixed 
potential  is  required).  Experimental  data  are  used  for  valida¬ 
tion.  Mass  transport  at  the  anode  is  a  factor  that  may  limit  the 
performance  of  the  fuel  cell. 

In  Sundmacher  et  al.  [46],  a  mathematical  model  was 
developed  to  describe  and  analyze  the  behavior  of  a  DMFC. 
The  main  phenomena  are:  convective  mass  transport  in  the 
anode  compartment;  methanol  oxidation  in  the  anode  cat¬ 
alyst  layer;  oxygen  reduction  in  the  cathode  catalyst  layer; 
methanol  and  CO2  mass  transfer  through  the  anode  diffusion 
layer;  methanol  crossover  and  undesired  methanol  oxidation 
at  the  cathode  catalyst  layer  as  a  result  of  methanol  crossover. 
A  multi-step  reaction  mechanism  (Eq.  (10),  in  which  the 
methanol  dissociative  chemisorption  was  assumed  to  be  the 
slow-step)  was  used  to  describe  the  electrochemical  methanol 
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Fig.  10.  Cell  voltage  response  to  periodically  pulsed  methanol  feed  con¬ 
centrations:  experimental  data  (top)  and  simulated  data  (bottom)  (reprinted 
from  [46],  Copyright  2001,  with  permission  from  Elsevier  Science). 


oxidation  at  the  Pt-Ru  catalyst: 

3Pt  +  CH3OH  Pt3  -  COH  +  3H+  +  3e“ 

3Ru  +  3H20  3Ru  -  OH  +  3H+  +  3e" 

Pt3  -  COH+2Ru  -  OH  Pt  —  COOH  +  H20+2Pt+2Ru 
Pt  —  COOH  T  Ru  —  OH  CO2  T  H2O  Pt  Ru 

(10) 


The  mixtures  at  the  anode  were  treated  as  a  pure  liquid 
phase  mixture.  To  complete  the  model,  mass  balances  for 
the  anode  compartment  and  the  anode  catalyst  layer,  and 
charge  balances  for  the  anode  and  cathode  catalyst  layers 
were  formulated.  The  model  of  the  effect  of  methanol  oxida¬ 
tion  at  the  cathode  considers  a  parasitic  current  proportional 
to  the  methanol  transport  through  the  membrane  (instanta¬ 
neous  methanol  oxidation  at  the  cathode).  Experimentally, 
the  authors  observed  that  the  cell  voltage  response  to  dynamic 
feeding  showed  a  significant  voltage  increase  after  a  sud¬ 
den  decrease  in  the  methanol  concentration  in  the  feed.  By 
simulation  studies  (Fig.  10),  it  was  possible  to  verify  that 
methanol  crossover  can  be  reduced  by  periodically  pulsing 
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the  methanol  feed  concentration,  resulting  in  an  increased 
cathode  potential  (and  an  increased  cell  voltage). 

The  authors  affirm  that  some  further  model  sophistications 
are  necessary,  such  as  considering  CO2  bubble  formation  in 
the  anode. 

In  Meyers  and  Newman  [47-49],  a  series  of  papers  regard¬ 
ing  the  simulation  of  a  DMFC  is  presented.  Modeling  of 
transport  and  kinetic  phenomena  is  considered.  The  model 
describes  the  multi-component  transport  of  species  in  the 
membrane  by  taking  into  account  a  thermodynamic  frame¬ 
work  developed  in  the  first  paper  of  the  series  [47] .  The  frame¬ 
work  provides  a  means  to  quantify,  in  the  second  paper  [48], 
the  gradients  in  electrochemical  potential  (for  hydrogen  ions, 
for  example,  /xH+  =  RgT  lnmH+  +  zH+  F<p)  for  species  in  the 
membrane  when  describing  the  driving  forces  for  the  multi- 
component  transport.  Besides,  a  kinetic  model  (considering 
reaction  steps  very  similar  to  those  presented  in  Nordlund  and 
Lindbergh  [42])  is  developed  to  describe  methanol  oxidation 
on  the  catalysts  [48].  In  the  third  paper  of  the  series  [49],  the 
mathematical  model  is  used  to  simulate  the  direct  methanol 
fuel  cell  and  quantify  general  aspects  of  its  design. 

In  Argyropoulos  et  al.  [50],  a  semi-empirical  model  based 
on  Tafel-like  kinetics  for  the  methanol  oxidation  and  oxygen 
reduction  reactions  and  on  effective  mass  transport  coeffi¬ 
cients  is  presented,  which  predicts  the  cell  potential  versus 
current  density  (Eq.  (11)): 

£cell  =  E%-  feeeii  log/  -  RI  +  C\  ln(l  -  C2I)  (1 1) 

In  Divisek  et  al.  [5 1  ] ,  a  two-dimensional  model  of  a  DMFC 
was  developed  and  tested.  The  model  considers  a  calculation 
domain  divided  into  different  zones,  with  the  diffusion  layers 
as  water-gas  systems.  The  authors  assume  that  the  transport 
of  water  follows  a  standard  two-phase  flow  mechanism.  This 
can  be  seen  as  a  generalized  Darcy’s  law.  The  pore  model 
for  the  transport  of  gases  describes  the  interactions  between 
the  gases  and  their  interactions  with  the  pore  walls,  incorpo¬ 
rated  into  the  two-phase  flow  modeling  approach.  Regarding 
the  dissolved  species,  a  standard  transport  model  is  consid¬ 
ered.  For  charged  species,  the  appropriate  potential  equations 
are  used.  Energy  considerations  are  derived  from  the  Fourier 
law,  with  convection  caused  by  the  flow  of  fluids  and  gases. 
The  reaction  models  account  for  the  methanol  (according  to 
Kauranen  et  al.  [52])  and  oxygen  kinetics.  Condensation  and 
evaporation  processes  are  considered.  The  final  equations  are 
discretized  in  time  and  space  with  the  objective  of  numeri¬ 
cally  solving  the  problem. 

In  Wang  and  Wang  [53],  two-phase  modeling  is  consid¬ 
ered  (such  as  in  Wang  et  al.  [22])  in  the  cathode  and  in 
the  anode  (due  to  CO2  evolution)  of  a  DMFC.  Tafel  kinetic 
equations  are  considered  to  describe  the  oxygen  reduction 
and  methanol  oxidation  reactions  (if  necessary,  adsorption 
steps  may  be  included  to  describe  methanol  oxidation).  The 
model  also  takes  into  account  the  mixed  potential  effect  of 
methanol  oxidation  at  the  cathode  as  a  result  of  methanol 
crossover.  Simulation  results  showed  (Fig.  11)  that  operat- 
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Fig.  11.  Polarization  curves  with  different  methanol  feed  concentrations 
(reprinted  from  [53],  Copyright  2003,  reproduced  by  permission  of  The 
Electrochemical  Society,  Inc.). 

ing  with  small  methanol  feed  concentrations  suffers  from 
low  limiting  current  densities.  In  the  medium  current  density 
range,  an  increase  in  methanol  feed  concentration  led  to  a 
small  decrease  in  cell  voltage  for  methanol  concentrations 
below  1 M.  For  a  feed  concentration  larger  than  2  M,  the  cell 
voltage  was  greatly  reduced  by  excessive  methanol  crossover 
and  the  maximum  current  density  was  limited  by  the  oxygen 
depletion  at  the  cathode  (resulting  from  excessive  parasitic 
oxygen  consumption). 

Additional  references  regarding  modeling  of  methanol 
fuel  cells  are:  Baxter  et  al.  [54]  (in  which  kinetic  and 
transport  limitations  are  studied);  Kulikovsky  [55]  (in  which 
some  local  features  due  to  geometric  aspects,  e.g.  the  low 
concentrations  of  reactants  in  front  of  current  collectors,  and 
mechanisms  of  methanol  transport  are  studied);  Fan  et  al. 
[56]  (in  which  some  local  features  generated  by  geometric 
aspects  are  also  studied);  Murgia  et  al.  [57]  (in  which  the 
effects  of  varying  some  empirical  parameters  are  examined). 
Some  works  consider  a  vapor  feed  DMFC,  such  as  Scott  et 
al.  [58]  (in  which  methanol  permeation  through  the  polymer 
membrane  and  the  caused  mixed  potential  are  studied); 
Dohle  et  al.  [59]  (in  which  the  impact  of  methanol  perme¬ 
ation  through  the  electrolyte  membrane  is  also  considered); 
Kulikovsky  et  al.  [60]  (in  which  a  new  type  of  current 
collector,  inside  the  backing  and  catalyst  layers,  which 
does  not  produce  shaded  regions  in  the  catalyst  layers,  is 
considered). 

Although  discussing  modeling  of  polymer  electrolyte 
membrane  fuel  cell  stacks  is  beyond  the  scope  of  this  paper, 
we  would  like  to  point  out  the  papers  of  Baschuk  and  Li  [61], 
Hamelin  et  al.  [62],  Jiang  and  Chu  [63],  Lee  and  Lalk  [64], 
Thirumalai  and  White  [65]  and  Amphlett  et  al.  [66].  Con¬ 
cerning  the  modeling  of  methanol  fuel  cell  stacks,  important 
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Fig.  12.  Reaction  diagram  of  the  microbial  fuel  cell  (reprinted  from  [73], 
Copyright  1995,  with  kind  permission  of  Springer  Science  and  Business 
Media). 

references  are  Argyropoulos  et  al.  [67-69],  Scott  et  al.  [70] 
and  Simoglou  et  al.  [71]. 

Considering  briefly  related  types  of  fuel  cells,  the  possi¬ 
bility  of  extracting  electricity  from  low  cost  and  abundant 
fuels  by  using  either  living  microorganisms  or  their  isolated 
enzymes,  in  biochemical  fuel  cells,  is  a  fascinating  point. 

In  a  microbial  fuel  cell,  the  direct  conversion  of  bio¬ 
chemical  energy  into  electricity  is  accomplished.  Basically, 
microbial  fuel  cells  are  divided  into  three  types  [72]:  the 
first  one  involves  the  direct  use  of  electroactive  metabolites, 
such  as  hydrogen,  produced  by  microbial  metabolism  from 
the  substrate;  the  second  type  involves  the  use  of  mediators 
for  electron  transport  from  a  certain  metabolic  pathway  to 
electrodes  (direct  contact  with  the  electrode  surface  is  not 
required);  the  third  one  involves  the  use  of  microorganisms 
that  may  grow  as  a  biofilm  on  the  electrode  surface  (media¬ 
tors  are  not  necessary).  When  only  enzymes  instead  of  living 
microorganisms  are  used,  in  vitro  systems  are  characterized. 

In  Zhang  and  Halme  [73],  a  bacterial  fuel  cell  consisting 
of  a  fuel  cell  device  and  a  bioreactor  is  considered  (Fig.  12). 

The  mediator  used  was  a  proper  quinone.  In  a  microbial 
fuel  cell,  besides  the  biochemical  processes  themselves,  there 
are  electrochemical  and  mass  transport  processes.  In  Zhang 
and  Halme  [73],  the  biochemical  processes  included  a  series 
of  enzymatic  reactions  that  could  be  modeled  by  Monod- 
like  equations.  Redox  reactions  were  considered  as  first  order 
processes.  All  mass  transport  processes  were  assumed  to  be 
fast  enough  when  compared  to  biochemical  and  redox  reac¬ 
tions.  The  process  could  be  described  by  considering  the 
proper  differential  equations  and  relationships,  such  as  the 
Tafel  equation  for  the  activation  overpotential.  The  simula¬ 
tion  studies  illustrated  how  the  current  output  depends  on 
the  main  system  variables.  The  effect  of  initial  substrate  con¬ 
centration,  for  example,  is  time  delayed.  A  suitable  mediator 
concentration  can  also  be  determined. 


3.  Conclusions 

Modeling  has  allowed  detailed  studies  concerning  the 
development  of  polymer  electrolyte  membrane  fuel  cells,  e.g. 
in  discussing  the  electrocatalysis  of  the  reactions  and  for  the 
design  of  water-management  schemes  to  cope  with  mem¬ 
brane  dehydration.  Results  show  that  for  cells  operating  with 
some  membranes,  water  transport  through  the  membrane  is 
a  limiting  effect;  for  others,  oxygen  diffusion  effects  (in  the 


gas  diffusion  electrode)  may  influence  the  potential-current 
relationship.  In  the  presence  of  hydrogen  dilution  in  the  fuel 
stream,  hydrogen  is  depleted  at  the  reaction  surface,  resulting 
in  substantial  anode  mass  transport  polarization  (which  can 
be  of  interest  when  considering  reformed  gas  as  the  anode 
feed).  Two-dimensional  models  have  been  used  to  represent 
reality,  but  three-dimensional  models  can  cope  with  some 
(important)  additional  aspects.  Consideration  of  two-phase 
transport  in  the  air  cathode  of  a  proton  exchange  membrane 
fuel  cell  seems  to  be  very  appropriate  (considering  a  non- 
isothermal  model  may  improve  the  analysis);  the  inclusion 
of  liquid  water  formation  and  transport  can  improve  the  pre¬ 
dictive  capability  of  the  model.  Knudsen  diffusion  has  been 
neglected  in  the  majority  of  the  papers. 

Regarding  direct  methanol  fuel  cells,  mass  transport  is  a 
factor  that  may  limit  the  performance  of  the  cell.  Adsorp¬ 
tion  theory  may  be  coupled  to  Tafel  kinetics  to  describe 
methanol  oxidation.  Methanol  crossover  must  also  be  taken 
into  account  (which  can  be  reduced  by  periodically  pulsing 
the  methanol  feed  concentration).  Extending  the  two-phase 
approach  to  the  DMFC  modeling  is  a  recent,  important  point. 
Two-phase  modeling  may  be  considered  both  in  the  cathode 
and  in  the  anode  of  a  DMFC. 

Regarding  other  types  of  fuel  cells  that  may  also  oper¬ 
ate  with  a  polymer  electrolyte,  the  possibility  of  extracting 
electricity  from  low  cost  and  abundant  fuels  by  using  living 
microorganisms,  in  biochemical  fuel  cells,  is  a  fascinating 
point.  In  contrast  to  the  increasing  experimental  approaches, 
few  authors  have  reported  simulation  studies  of  this  kind 
of  system.  The  study  of  biochemical  fuel  cells  by  applying 
proper  modeling  techniques  (coupled  to  experimental  stud¬ 
ies)  can  help  developments  in  this  area,  and  certainly  should 
be  more  explored  in  the  future. 
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